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-q. ABSTRACT 

O 
O 

■ Galaxy disk formation must incorporate the multiphase nature of the 



interstellar medium. The resulting two-phase structure is generated and 
Q ■ maintained by gravitational instability and supernova energy input, which 

yield a source of turbulent viscosity that is able to effectively compete in the 
protodisk phase with early angular momentum loss of the baryonic component 
via dynamical friction in the dark halo. Provided that star formation occurs on 
^ ■ the viscous drag time-scale, this mechanism provides a means of accounting for 

disk sizes and radial profiles. The star formation feedback is self-regulated by 
O ■ turbulent gas pressure-limited percolation of the supernova remnant-heated hot 

phase, but can run away in gas-rich protodisks to generate compact starbursts. 
O ■ A simple analytic model is derived for a Schmidt-like global star formation law 

^ ■ in terms of the cold gas volume density. 

S ' Subject headings: galaxies: format ion-galaxies: ISM-galaxies: kinematics and dynamics- 

■ galaxies: stellar content 

> 
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5-h ' 1. Introduction 

An outstanding problem in understanding the formation of galaxy disks is the origin 
of the characteristic disk scale length. In the absence of a theory that determines this 
parameter, one can place little credence in simulations of numerical or semi-analytical 
galaxy formation. Yet until recently, the elements of a robust theory for disk scales appeared 
to be in place. 

Tidal torques between neighbouring density fluctuations generate an initial 
dimensionless angular momentum Aj = J\E\G~ l M~ 5 / 2 rs 0.05, with a broad dispersion 
AAj ~ Aj. Baryons cool and contract once dark halos virialize, conserving specific angular 
momentum to form a centrifugally supported disk when A ~ 1. The final disk radius is 
r d ~ \ r hi where is the halo radius. This yields a disk scale length of several kpc for an 
galaxy, as observed (Fall and Efstathiou 1980). 
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However with the advent of numerical hydrodynamical simulations with adequate 
resolution, it has become apparent that angular momentum is effectively transferred 
outwards as dense baryon clumps sink into the central regions of the forming galaxy, 
dumpiness in the dark halo also tidally perturbs the disk and aggravates the angular 
momentum transfer problem. The result of the simulations is that the disk sizes are 
invariably too small (Steinmetz and Navarro 1999), by a factor ~ 5. 

Two schemes have been proposed in order to halt excessive angular momentum transfer. 
Photoionization of the intergalactic medium destroys small-scale substructure (Navarro 
and Steinmetz 1997). These structures are the building blocks of large-scale structure and 
a source of dumpiness during massive halo formation. However the CDM structure is 
unaffected and there is little effect on the formation of massive galaxies. Another problem 
related to that of angular momentum transfer in that it involves baryonic dissipation is that 
of overcooling. This results in premature exhaustion of the gas supply. A solution designed 
to tackle both of these problems appeals to feedback from star formation via supernovae. 
Provided gas cooling is delayed to z ~ 1, one can form disks of the appropriate size and 
mass (Weil, Eke and Efstathiou 1998). 

However this success comes at a price. Late disk formation implies strong evolution in 
disk size, luminosity and surface brightness at z ~ 1. Large disks are observed to z ~ 1, 
and there is little evidence for any evolution in scale length and surface brightness, other 
than passive stellar population ageing (Lilly et al. 1998). Disks are observed to z ~ 1 with 
a Tully-Fisher relation that is essentially indistinguishable from that of present-day disks, 
apart from passive evolution. Not all analyses agree: Bouwens and Silk (2000) find that 
up to 2 mag of surface brightness evolution is required at z ~ 1 in an analysis of a much 
larger data set (Simard et al. 1999). However the very existence of substantial numbers of 
disk galaxies at z ~ 1 means that early heating of protodisks and consequent suppression of 
star formation to z ~ 1 appears to be an unacceptable resolution of the angular momentum 
transfer problem. 

A new approach is required. I have previously proposed (Silk 1997) the following 
model for disk formation. Feedback from supernovae results in a two-phase interstellar 
medium, the porosity of which controls feedback on the supply and compression of cold 
gas and the ensuing star formation rate. Gas-rich disk gravitational instabilities drive 
cloud-cloud collisions that provide an effective turbulent viscosity which is responsible for 
the formation and contraction of the disk. Viscous galactic disk models have a long history, 
commencing with the original proposal by Silk and Norman (1981), the definitive study 
by Lin and Pringle (1987a) that derived exponential stellar surface density profiles, and 
the chemical evolution studies of Clarke (1989,1991). In the cosmological context of cold 
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dark matter halos, turbulent viscosity has been incorporated into disk models by Firmani 
et al. (1996). I show here that the viscous drag time-scale is of the same order as the 
dynamical friction time-scale that controls the loss of angular momentum from baryons to 
dark matter (Section 2). I apply porosity feedback to derive a Schmidt-like star formation 
law that depends on cold gas volume density (Section 3), and I examine outflows and radial 
flows (Section 4). The requirement that the star formation time-scale is on the order of 
the viscous time-scale suffices to give an exponential disk with an appropriately large-scale 
length and that moreover forms at a relatively early stage (Section 5). 

2. Angular momentum problem resolved 

The loss of baryonic angular momentum during disk formation is primarily due to the 
dynamical friction exerted by the inner dark halo on the baryonic clumps. The time-scale 
for dynamical friction to act is approximately (M / M c i OU( i) dynamical times, or 

«*-(f)(s)"-fO" 

where the first equality made use of the tidal radius l t defined by 

p/p = {k/lf 

for a cloud of density p and size / in a medium of mean density p. The second equality 
replaces the disk circular velocity V by fir and a g by Ql t , the latter being the cloud velocity 
dispersion induced by gravitational instabilities (Gammie, Ostriker and Jog 1991). 

To compare with the viscous drag time-scale, defined by 

t v = r 2 / v, 

I compare three alternate definitions of viscosity. In a non-self-gravitating cold accretion 
disk, the a-prescription (Shakura and Sunyaev 1973) is commonly used: 

v\ = aHc s a ~ 1, 

where H is the scale height and c s is the sound velocity. This prescription is lacking 
in fundamental motivation. A recent proposal (Duschl, Strittmatter and Biermann 
2000) appeals to turbulence experiments which suggest that a limiting Reynolds number 
R cr it ~ 10 2 — 10 3 demarcates the onset of turbulence, and yields the ansatz 



v 2 = (3rV : (3 « R~ 
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This expression should be applicable to self-gravitating systems and to spheroidal as well as 
disk geometries. Finally, non-axisymmetric gravitational instabiliies of a cold disk provide 
an effective turbulent viscosity that transfers angular momentum (Lin and Pringle 1987b). 
This results in a characteristic instability scale (It) and turbulent velocity (a g ~ Qlt), 
motivating the hypothesis that 

z/ 3 = 'ykcTg, 7 ~ I. 
The corresponding time-scales for viscous drag are 

In all three prescriptions for viscosity, the viscous time is of order the dynamical friction 
time; and indeed t vz (and even t U2 ) can be substantially less. It is likely that only t U2 and t U3 
are relevant to self-gravitating disk formation, and I conclude that viscosity in the protodisk 
can operate more effectively than dynamical friction. 

This conclusion would not resolve, and indeed would even aggravate, the angular 
momentum problem, since viscous drag also removes angular momentum from the gas 
clouds, were it not for the fact that the star formation time t* itself likely to be comparable 
to the viscous time. This is plausible because the radial gas flow M g = 27rS 3 z/ is regulated 
by viscosity, and the gas flow controls the gas supply available for star formation at any 
given radius. Strong support for this conjecture comes from the model of a self-gravitating 
viscous disk, in which Lin and Pringle (1987a) find that t v results in an exponential 
surface density profile. Hence I conclude that star formation develops before substantial 
angular momentum loss is induced by dynamical friction of infalling clouds against the 
dark matter, and that the resulting disk scale- length is therefore of order XiTh, where 
is the initial halo radius (~ 100 kpc) at overdensity 200 (the virialization radius) and 
Aj ~ 0.05 is the dimensionless halo angular momentum parameter induced by gravitational 
torques in the early universe. Moreover the resulting stellar disk surface density profile is 
exponential. Of course, numerical simulations are needed to verify these conjectures and 
should eventually be feasible. 



3. Global star formation rate: the second parameter 

It is usually assumed that the star formation rate in a self-gravitating gas-rich disk is 
regulated by gravitational instability according to whether the Toomre parameter 

Sna„ 
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is smaller or larger than unity, for instability or stability respectively. Here, k is the 
epicyclic frequency, a g is the gas-velocity dispersion, S 9 is the gas surface density, and 5 is 
a parameter, of order unity for a gas-rich disk, that allows for the stellar contribution to 
the self-gravity of the gas component of the disk. Ha, HI and CO mapping of disks shows 
that Q g ~ 1 in the star-forming regions of disk galaxies, suggesting that star formation is 
self- regulated by gravitational instability (Kennicutt 1989). This is plausible, since disk 
instabilities heat the disk and raise Q g via tidal acceleration of clouds, whereas gas cloud 
collisional dissipation decreases the gas velocity dispersion, decreasing Q g . Star formation 
occurs in the clouds that form by gravitational instability, grow in mass by coalescence, and 
become unstable to fragmentation. Q g self-regulation predicts that there is a threshold to 
star formation. Star formation ceases when the gas surface density drops below a critical 
value 

Ecrit ~ 5K<J g /7TG. 

Now k ~ 2^f2 for a flat rotation curve, f2 « 10~ 15 s _1 and a g ~ 10 kms -1 for the Milky 
Way, as a typical example, so that T, crit ~ 10M o pc -2 . Clouds of giant HII regions, both in 
our galaxy and in nearby spirals, support this star formation threshold (Kennicutt 1998). 
However Ha mapping (Ferguson et al. 1998) has shown that substantial star formation 
occurs at E < E crit . To properly account for star formation, a second parameter is evidently 
needed that supplements or modifies the Q g criterion. 

A promising direction for such a modification would involve the dependence of the star 
formation rate on gas volume density, rather than surface density. Disk flaring develops in 
the outer disk, and the resulting density decrease at fixed E 5 would be a logical inhibitor of 
star formation that occurs near, but independently of, the radius at which E g drops below 
E crit . Indeed, it is entirely possible that the physics that controls disk flaring, namely the 
transition between the dominance of disk and halo self-gravity, is ultimately responsible for 
determining the transition between the Q g < 1 (inner disk) and Q g > 1 (outer disk) regimes. 
I develop here a physical criterion for controlling the star formation rate via feedback that 
depends primarily on the gas density, rather than the gas surface density. Supernovae 
control the energy and momentum input that heats and stirs the interstellar medium. The 
rate of supernovae is controlled by the star formation rate, and therefore provides a natural 
means of self-regulating star formation. The essence of this self-regulation can be described 
by a two-phase model of the interstellar medium, consisting of the hot gas directly heated 
by supernova remnants and the cold gas that is compressed by supernova blast waves into 
dense shells, dense cold clumps of which also survive the passage of the supernova blast 
waves. 

I have argued that global feedback operates via the porosity of a two-phase interstellar 
medium (Silk 1997). I now further extend these arguments by deriving an expression for 
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the star formation rate that does not explicitly depend on Q. The porosity is related to the 
hot gas fraction by P = — ln(l — /), where / is the volume filling factor of hot gas. For 
global feedback, the hot gas volume fraction must be significant, say / > 0.1. If the gas 
density or pressure were sufficiently high, supernova feedback would be restricted to such 
small volumes that there would be no global effect. A two-phase interstellar medium is 
adopted in order to describe the interaction of cold and hot gas. The cold gas is atomic and 
molecular gas, and initially comprises the entire disk. Once stars form and die, a steady 
state is reached in which bubbles of hot gas driven by supernovae sweep out and envelop 
much of the cold gas. The hot gas mass fraction is determined by cold cloud evaporation 
and by cooling of the hot gas. The diffuse hot gas erodes cold clouds and leaks out of the 
disk via chimneys and fountains. I will estimate below the rate at which gas ejection occurs. 
The porosity is a convenient and calculable parameter that defines the filling factor of the 
hot phase, and hence its ability to diminish the efectiveness of star formation by eroding 
the cold gas supply. 

Supernova remnants deposit energy and, ultimately, momentum into the interstellar 
gas at a rate p\v S N where vsn is the specific momentum injected by supernovae per unit 
star formation rate p\ and is given by 



conserving regime, E 51 = _E S 7v/10 51 ergs the supernova energy (taken to be 10 51 ergs), and 
( g the metallicity relative to solar of the ambient gas (Cioffl, Mckee and Bertschinger 1988). 
Also msN is the mean mass required in forming stars in order to produce a supernova. For 
SNII, one simply assumes an initial mass function (IMF) with all stars of mass above 8 
M becoming supernovae, so that tusn ~ 250 M for a Miller-Scalo IMF. One can crudely 
double the inferred rate for an estimate of the rate of Types I and II supernovae after the 
first 10 8 years have elapsed (to allow sufficient time for SNIa to form). 

Porosity is defined to be the product of the supernova remnant 4-volume at maximum 
extent, when halted by ambient gas pressure, and the rate of bubble production. The 
porosity is a measure of the fraction of volume / occupied by the hot phase (T ~ 10 6 K) 
associated with the interiors of supernova remnants. By definition f = 1 — e~ p , and the 



stars formed per supernova, and R a is the radius of the supernova remnant at time t a when 
halted by the ambient (turbulent) gas pressure p g . One finds that the porosity 



v S n = E S N/v c m S N = 500Ell /14 n g 1/7 m 2 5oCg 3/14 kms 1 



(1) 



with v c = 413£j[ 4 nJ/ 7 C* /14 b 





where p\ is the star formation rate, msN is the mass in 



P oc p* P, 



is extremely sensitive to the interstellar pressure. This provides the motivation for the 
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feedback prescription. By incorporating an analytic fit to the evolution of a spherically 
symmetric supernova-driven shell (Cioffl et al. 1988), one can write 

P = G-\af V 3 V' n P* , (2) 

where p* is the star formation rate, p g is the ambient gas pressure, both thermal and 
turbulent, and aj is a fiducial velocity dispersion that is proportional to E^m^(~ - 2 and 
may be taken to be 22 km s _1 for E$n — 10 51 erg, tusn = 250 M Q and ( g = 1. Note that at 
large pressure, ?< 1, and porosity is primarily controlled by the ambient pressure, which 
I take to be dominated by turbulence: p g = p g a 2 . 

Feedback is regulated by the porosity parameter P. Previously, the star formation rate 
had no explicit feedback. Consider the possibility of strong feedback. In this case, the filling 
factor of hot gas is of order fifty percent, which is the requirement for strong feedback. 
Rewriting (2) as 

P^G^pf (^fV (3) 

one explicitly incorporates feedback into the star formation rate. The natural limit in which 
feedback self-regulates has P ~ / ~ 1/2. If P <C 1, feedback is irrelevant: if P 3> 1, star 
formation would be completely suppressed. It seems likely, though hard to prove rigorously, 
that / ~ 1/2 corresponds to self-regulation of disk star formation. Certainly, / ~ 0.2 — 0.5 
prevails in the local interstellar medium (Slavin and Cox 1993; Ferriere 1995). In this case, 
the star formation rate is independent of Q, and is controlled by ambient pressure. The 
turbulent pressure is fed by disk gravitational instabilities, and implicitly depends on Q. 
The derived global star formation law has a normalisation, equivalent to star formation 
efficiency, that depends explicitly on cloud velocity dispersion. The star formation rate and 
efficiency are enhanced by stirring: this might be expected near a bar or as a consequence 
of a merger. 

Another limiting case is P ~ / << 1. Here porosity is small and feedback is 
unimportant. This might be relevant in a dense gaseous protodisk, and would suggest the 
possibility of runaway star formation. 



4. Feedback and outflow 

Supernovae are responsible for the feedback from star formation that renders star 
formation inefficient. The star formation efficiency e may be defined explicitly by writing 
the star formation rate as 

p\ = eflp g , 
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where e is the fraction of gas converted into star per dynamical time. Adapting the porosity 
model for feedback, one now can write the porosity as 




where p is the total (gas plus star) density in the disk or protodisk. High star formation 
efficiency drives high porosity and significant feedback. 

The momentum input from supernovae is dissipated via cloud-cloud collisions and 
outflow from the disk. In a steady state, the momentum input rate p*vsn must balance 
the cloud collisional dissipational rate pgcrpy 1 and the momentum carried out in outflows 
fPgH' 1 , where p g is the turbulent pressure p g <jg of the two-phase interstellar medium, 
H = <jg/27rGY> is the disk gas scale height, and £ is the surface mass density. At a given 
disk radus, there is also viscosity-driven inflow, M g = 27rS 5 ^, which provides a gas supply 
for star formation. Ignoring the outflow and infall contributions in the momentum budget, 
one has 

eVtpgVsN = PgCTgQ, 

so that e = 0.02 (cr g /10 kms -1 ) (500 kms -1 /vsn)- This suggests that supernova feedback 
can indeed yield the required low efficiency of star formation. 

For a galaxy such as the Milky Way, the global star formation efficiency is expected to 
be around 2 percent, both as inferred from the global values of gas mass (~ 6 x 10 9 M Q ) and 
star formation rate (~ 3M q jt~~ 1 ) after allowance for gas return from evolving stars (the 
returned fraction ~ 0.5 for a Miller-Scalo IMF) over a galactic dynamical time and as more 
directly inferred from studies of HII region radio luminosities summed over molecular cloud 
masses (e.g. Williams and McKee 1997). 

Inserting the derived expression for star formation efficiency into the equation for the 
porosity, I find that 

P „ feV' 7 in. = o JeM^^^y" . 

W W v sn V 0- 1 / V10kms-V 

The observed three-dimensional cloud velocity dispersion is 11 kms -1 (for molecular clouds 
within 3 kpc of the sun) (Stark and Brand 1989). Thus the Milky Way interstellar medium 
has predicted porosity comparable to what is observed. 

If the hot volume fraction is large, outflow from the disk is inevitable. The ratio of 
outflow to star formation rate per unit area of disk be written as (f s p g H~ 1 V~^)(eQp g )~ 1 , 
where V esc is the escape velocity from the disk. Now p g = 7iGT, g a g = ^H g QV and 
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Tj g = 2p g H, so that the outflow rate simplifies to 

outflow rate f s 
star formation rate e ' 

where f s is the surface covered in hot bubbles. 

This ratio will typically be, for say / ~ P ~ 0.3 and e ~ 0.03, around 10, and is 
comparable to what is observed for both star-forming dwarfs and more luminous disk 
galaxies (Martin 1999). In general, the bulk of the gas cannot escape in a wind from normal 
galaxies, otherwise the gas reservoir would be seriously depleted. Known infall rates amount 
to at most 10 percent of the star formation rates. The gas must cool in the halo and fall 
back into the disk. 

The viscosity prescription leads to a radial inflow of the gas. This amounts to 
v r = z/r" 1 , which can be rewritten as 

v r = jaghr' 1 = JCTg/V, 

where I have used the expression for a g derived for gravitational instability-driven random 
motions. Hence one typically expects 

v r « 0.57 ( — ° 9 _\ ( s — jkms" 1 . 



.lOkms -1 / V V 
The associated mass flux is 

M g = 2iruT,g = 27T7/i<7 g E 5 . 
Again making use of a g = flit and £ = QV/2ttG, the mass flux reduces to 

M g = t^VG" 1 (£,/£) = 7 K/l0 2 fiM(£,/£). 

Comparing this with the global star formation rate M* = eQM g , with e = <j 9 /vsn, one 
obtains 

M g = 1 {<j g /V){v SN /V)M*, 

or 

m 9 /m* = (1) ( ggg r ) faooto^y 

91 \8J VlOkms" 1 / V500kms~V \ V J 

I conclude that the viscous supply of gas may account for of order 30 percent of the net gas 
supply, after return, as is required to balance star formation. The associated radial flow is 
sufficient to modify chemical evolution in the solar neighbourhood and resolve the G dwarf 
problem (Clarke 1991), and to generate a disk metallicity gradient (Clarke 1989). However 
the most novel application is to the disk scale length, as I now demonstrate. 



-10- 



5. Disk properties 

The disk size may be estimated by setting the viscous time-scale equal to the star 
formation timescale. This guarantees an exponential profile, and the characteristic length 
scale is obtained by setting t v = t s f r , so that 

Td,u = J(Tgt s f r /V. 

Inserting the star formation time-scale 

eQ Pg ft \pj \a g J 

I infer that 

_ f lVsNcr g \ (p_ \ _ f v,sN \ ( o g \ ( 5OM pc" 2 \ / _p_\ 
^ - \ 2vrGE ){pj ~ Z1 U00kms-V vlOkms-V ^ E ) \2pJ PC ' 

comparable to the scale lengths of Milky Way-type disks. 

It is interesting to compare this with the standard derivation of the final radius of a 
nonlinear cold dark matter density perturbation that collapses to approximate virialization 
at a spherically-averaged overdensity of 200. The radius at mean overdensity 200 is 

0.011/ 

r 200 ~ F) 

where V is again the circular velocity, H(z) is the Hubble parameter, and fto is the density 
parameter. The baryons dissipate and contract within the cold dark matter halo to a final 
disk radius, if angular momentum is conserved, of 

r d,cdm — \ r 2Q0 

In the viscous model, the disk radius can be written in terms of the virial radius defined by 
the self-gravitating baryonic component, r v = V 2 (ttGT i )~ 1 , so that 

_7 ( VSN(T 9 \ P 

d ' v 2 \ V* J p 9 v - 

Allowing for the initial baryon fraction ft, (expected to be ftb/ft, m ~ 0.05/0.3 ~ 1/6), one 
finds that the gas first becomes self-gravitating at radius ^200- Hence the disk radius is 
comparable to that attained in the idealized cold dark model, with conservation of specific 
angular momentum of the baryons, 



rd,u _ fbl ( V SNVg \ ( P_ N 

r d ,cdm \ 2 V V 2 J \p 9/ 



1. 
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The strong dependence of r v (or r 20 o) on z leads to the predicted sensitivity of 
the Tully-Fisher relation to redshift (Mao, Mo and White 1998), in contradiction with 
observations (Voigt 1999). In the viscous model, this effect is diluted because of the 
proportionality of disk radius to gas velocity dispersion. I expect a g to increase with 
redshift because the effect of mergers and infall will result in systematically more turbulent 
protodisks. Hence this should partially compensate the reduction of at earlier epochs 
arising from the explicit redshift dependence in the expression for r 20 o- Some reduction inr^ 
is neded in order to account for the observed evolution in disk surface brightness at z ~ 1 
(Bouwens and Silk 2000). 

The Toomre parameter for disk instability may be rewritten as 




where £ is the total disk surface density of stars plus gas. Initially Q g < 1, so that the 
protodisk was unstable. At late times however, Q g increases as £ 3 decreases, although 
o g also decreases, so this increase in Q g is partially compensated by the decrease in gas 
turbulence, and the disk should remain at least marginally unstable, Q g < 1. 

To maintain Q g ~ 1 over many dynamical times, as required for disks of spiral galaxies, 
the disk must remain gas-rich. Gas cooling is required to dissipate the dynamical instability 
heating. One can estimate that the cold gas controls the gravitational instability of the 
disk if H g /cr g > /L^/cr*, and this condition requires a gas fraction of 30 percent or more, if 
the gas velocity dispersion a g ~ 10 km s _1 and the stellar velocity dispersion ~ 30 km s -1 . 
This is indeed observed for the solar neighbourhood, where fi g ~ 15 M pc -2 and //* ~ 40 
M pc -2 . The gas supply is prescribed by radial inflow in the viscous disk model, as is 
consistent with the gas reservoir inferred from observations of the high velocity clouds. 



6. Discussion 

In summary, I have argued that the justification for the viscous disk model comes from 
the fact that viscous and star formation time-scales are likely to be comparable, 

hlL = ry (JL) (°a) ( VsN \ 

t» ' \Pg) \Vj V1W' 

where Qd is the rotation rate evaluated at the disk scale length. The baryons retain most of 
their initial angular momentum since the dynamical friction time-scale exceeds the viscous 
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time-scale, 

£ = •*(-)■ 

This consequently justifies an exponential disk profile with scale length of order 
r d,u ~ fb(vsNC g /V 2 )r2oo- This scale coincidentally happens to be of order A^ooj confirming 
that approximate angular momentum conservation is achieved, and disk sizes match the 
observed range. 

The viscous disk model predicts a radial inflow rate of order 10% of the disk star 
formation rate. This will generate a gradient in the chemical abundances, and also provides 
a source of gas for bulge formation during the phase of early gas-rich disk evolution. The 
predicted ratio of bulge to disk stellar mass is of order cr g vsN/V 2 , i-e. of order 10%. Late 
forming bulges by early viscous disk evolution would imply that bulges are more metal-rich 
than the disk, and that the age spread of bulge stars is comparable to that of the old disk. 

Supernova feedback regulates the star formation rate. The dominant parameter is 
porosity, largely determined by the turbulent gas pressure p g . When p g is in the range 
expected for typical disk galaxies, P ~ 1 and feedback is important. I assume that the 
feedback is negative: further work is needed to justify this. In the case of high turbulent 
pressure, expected in the aftermath of mergers, Pel. The lack of feedback means that 
star formation can run away, only being limited by the available gas reservoir. In general, I 
find that supernova feedback results in a star formation rate 

» ttp g (a g /v SN ). 

Thus even in starbursts, a Schmidt-type relation is maintained. The predicted efficiency 
can be as high as 20 %, for merger-induced turbulence a g ~ 100 kms -1 and a standard 
IMF, for which vsn ~ 500 kms" 1 . 

A top-heavy IMF could substantially reduce vsn, and star formation efficiences of 
order 50% would then be attainable. Such efficiencies may be needed in order to account 
for the luminosities measured in many ultraluminous infrared galaxies where the molecular 
gas masses are measured. A top-heavy stellar initial mass function might be required in 
protogalaxy mergers in order to reconcile the hypothesis that ellipticals form in such events 
with the observed paucity of young ellipticals at intermediate redshifts. 

Gas density (and pressure) rather than surface density is the dominant gas parameter 
in porosity-regulated star formation. This means that disk flaring at the transition from 
disk self- gravity to halo gravity dominance, where the gas scale height increases due to the 
reduction in disk self-gravity, will tend to result in a reduction of the disk star formation 
rate. Disk gravitational instability, which controls the supply of massive molecular clouds, is 
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also quenched at approximately the same radius. This leads to a more complex dependence 
of star formation rate on gas surface density than suggested by the simple Q g threshold 
models that have hitherto been advocated. 

I thank J. Devriendt for discussions and a critical reading of the manuscript. While 
this paper was in a refereeing phase that was protracted for reasons beyond my control, a 
paper appeared that contains an extensive discussion of feedback in a similar context to 
that discussed in the present paper (Efstathiou 2000). A more detailed comparison of the 
two approaches will be presented elsewhere. 
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